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ABSTRACT 

We argue that the observed correlation between the radio luminosity and the thermal X-ray luminosity 
of radio emitting galaxy clusters implies that the radio emission is due to secondary electrons that 
are produced by p-p interactions and lose their energy by emitting synchrotron radiation in a strong 
magnetic field, B > (SnaT^^jg) 1 ^ 2 ~ 3 /xG. We construct a simple model that naturally explains the 
correlation, and show that the observations provide stringent constraints on cluster magnetic fields 
and cosmic rays (CRs): Within the cores of clusters, the ratio (3 COIC between the CR energy (per 
logarithmic particle energy interval) and the thermal energy is /3 coro ~ 2 • 1CP 4 ; The source of these 
CRs is most likely the cluster accretion shock, which is inferred to deposit in CRs ~ 0.1 of the thermal 
energy it generates; The diffusion time of 100 GeV CRs over scales > 100 kpc is not short compared 
to the Hubble time; Cluster magnetic fields are enhanced by mergers to > 1% of equipartition, and 
decay (to < 1 fiG) on 1 Gyr time scales. The inferred value of /3 COT e implies that high energy gamma- 
ray emission from secondaries at cluster cores will be difficult to detect with existing and planned 
instruments. 

Subject headings: acceleration of particles - galaxies: clusters: general - radiation mechanisms: non- 
thermal - radio continuum: general - X-rays: general 



1. INTRODUCTION 

Nonthermal emission is observed in several clusters of 
galaxies, mainly in the radio band: giant radio halos 
(RHs) and mini-radio halos, fairly symmetric sources 
at the cluster center, as well as radio relics, elongated 
sources at the cluster periphery, have been observed (e.g., 
Feretti & Giovannini 2008). The radio emission is inter- 
preted as synchrotron radiation, thereby suggesting that 
relativistic electrons and magnetic fields are present in 
the intra-cluster medium (ICM). While the radio relics' 
emission can be attributed to merger or accretion shock 
waves (e.g., Ensslin et al. 1998), the origin of the radio 
halos is not yet understood. 

Pointed radio observations of a complete sample of 
about 50 X-ray clusters at rcdshift z = 0.2 — 0.4, 
with rest-frame [0.1,2.4] keV luminosity Lx[0. 1,2.4] > 
5 • 10 44 ergs _1 , have been recently carried out (Venturi 
et al. 2007; Brunetti et al. 2007; Venturi et al. 2008). 
One of the findings of these observations is that ~ 30% 
of the X-ray luminous clusters host RHs, with all RHs 
being in merging clusters. In addition, a tight correlation 
between the radio luminosity and the X-ray luminosity 
of clusters with RHs was established. Note, that only a 
few halos have good multi-frequency observations, that 
allow one to estimate their integrated spectrum. These 
spectra are usually consistent with a spectral index of 1, 
L v cx v^ 1 (Feretti & Giovannini 2008). 

Several models for the synchrotron emission in RHs 
have been presented in the literature. These models dif- 
fer in the assumptions regarding the origin of the emit- 
ting electrons. In some models the emitting electrons are 
secondary electrons and positrons that were generated by 
p-p interactions of a CR proton population with the ICM 
(e.g. Dennison 1980; Blasi & Colafrancesco 1999) while in 
others the emitting electrons are reacelerated by turbu- 
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lence from a preexisting population of nonthermal seeds 
in the ICM (secondary or otherwise, e.g. Brunetti et al. 
2001; Petrosian 2001; Brunetti & Blasi 2005; Cassano & 
Brunetti 2005; Cassano et al. 2007; Brunetti et al. 2008). 

In this paper we present a simple analytic model, fol- 
lowing Kushnir & Waxman (2009b), for the radio emis- 
sion in clusters with RHs, and show that it naturally 
reproduces the observed correlation between the radio lu- 
minosity and the thermal X-ray luminosity in such clus- 
ters. The observations are described in § 2, and the sim- 
ple model for the X-ray and radio emission from clusters 
with RHs is described in § 3. In this model, the radio 
emission is produced by secondary electrons, that are 
produced via p-p interactions of intra-cluster CRs with 
the ICM, and lose their energy by emitting synchrotron 
radiation in a strong magnetic field. It is assumed that 
the magnetic field is strong enough to ensure that syn- 
chrotron losses dominate the electrons' energy loss, and 
that the loss time is short compare to the cluster dynam- 
ical time. The validity of these assumptions as well as 
the implications of the bimodality of the radio luminosity 
distribution are discussed in § 4. In § 5 we show that al- 
ternative models for the radio emission, with lower values 
of the magnetic field or different sources for the emitting 
electrons, are unlikely to reproduce the observed corre- 
lation between X-ray and radio luminosity. Our results 
are summarized and their implications are discussed in 
§6. 

We note that in some clusters nonthermal emission 
is also observed in hard (> 20keV) X-rays (HXR) (for 
review, see Rephaeli et al. 2008). The HXR emission 
is usually interpreted as due to inverse Compton scat- 
tering of cosmic microwave background (CMB) photons 
by nonthermal relativistic electrons (e.g., Rephaeli 1979; 
Sarazin 1999). We discuss the origin of the nonthermal 
HXR emission in a separate paper (Kushnir & Waxman 
2009a). We show there that both the HXR observations 
and the results of the analysis presented in this paper, 
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Fig. 1. — The distribution of clusters in the Pi. 4 — Axfo.1,2.4] 
plane, taken from Brunetti et al. (2007) (reproduced by permission 
of the A AS). Over plotted (black dashed line) is the best linear fit, 
eq. (1), for clusters that contain RHs (see text). 



which indicates that /3 COIe ~ 2 • 1CP 4 , imply that the 
nonthermal HXR emission can not be produced by the 
same electrons producing the RHs. Both HXR and RH 
emission are consistent, instead, with a model where the 
nonthermal particles responsible for both radio and X- 
ray emission originate in cluster accretion shocks. In this 
model, the HXR emission is due to IC energy loss of elec- 
trons accelerated directly by the accretion shocks, while 
the radio emission is due to synchrotron emission of sec- 
ondary electrons produced by the interaction of cosmic 
ray (CR) protons, which are also accelerated in the accre- 
tion shocks, with the thermal ICM protons (see Kushnir 
& Waxman 2009b, a, for detailed discussion). 

2. OBSERVATIONS 

We focus on the distribution of the clusters in the ra- 
dio and X-ray luminosity plain. Pointed radio obser- 
vations of a complete sample of about 50 X-ray lumi- 
nous (L x [0. 1,2.4] > 5 • 10 44 ergs _1 ) clusters at redshift 
z = 0.2 — 0.4 have been recently carried out (Venturi 
et al. 2007; Brunetti et al. 2007; Venturi et al. 2008). 
The distribution of clusters in the Pi. 4 — L X [o. 1.2.4] plane 
is shown in figure 4 of Brunetti et al. (2007), where Pi. 4 
is the radio luminosity at 1.4 GHz. This figure is repro- 
duced in figure 1. 

The main characteristics of the distribution can be 
summarized as follows: 

1. The radio luminosities of most of the clusters 
that contain RHs follow a linear correlation be- 
tween ln(Pi.4) and ln(Lx[o.i,2.4])- The disper- 
sion around the correlation for these clusters 
is relatively small: The standard deviation of 

ln(Pi.4,observed/Pi.4,e X pected) is ~ ln(1.7). Since this 
dispersion is much larger than the measurement er- 
rors, it reflects an intrinsic scatter of the radio lu- 
minosity (around the linear correlation). 

2. The best linear fit, obtained assuming that 
the intrinsic variance around the correlation is 



Aln(Pi. 4 ) = ln(1.7), is 2 

P1.4 * 5 • lO^xfo.^.^ergs- 1 Hz' 1 , (1) 

where £x[o.l,2.4],45 = £x[o.i,2.4]/10 45 ergs" 1 . 

3. About ~ 1/3 of the clusters contain RHs. For 
bright clusters, Lx\o. 1,2.4] ^ 10 ergs -1 , that do 
not contain a RH, the upper limit on Pi. 4 is about 
an order of magnitude smaller than the value im- 
plied by the correlation. 

For a few clusters, we verified that the RH luminosi- 
ties represent the total radio luminosities of the clusters 
by studying their radio surface brightness radial profiles 
(taken from Cassano et al. 2007; Govoni et al. 2001). 

As we show below, the observations described above 
allow us to derive stringent constraints on the magnetic 
fields and on the CR population within clusters. 

3. A SIMPLE MODEL FOR RADIO EMITTING 
CLUSTERS 

In this section we present a simple model for the non- 
thermal emission of clusters that have RHs, and show 
that it naturally reproduces the observed correlation be- 
tween the radio and X-ray luminosities. In this model 
intra-cluster CRs interact with the thermal gas of the 
ICM to produce secondary electrons and positrons via 
p-p interactions. We assume that secondary electrons 
and positrons lose their energy primarily by emitting syn- 
chrotron radiation in a strong magnetic field, B, and that 
the energy loss time is short compared to timescales over 
which the magnetic field or the CR electron injection 
change considerably. Under these assumptions, the va- 
lidity of which is discussed in § 4, the radio synchrotron 
luminosity is equal to the energy generation rate of sec- 
ondaries through p-p interactions. This generation rate 
is in turn correlated with the thermal Bremsstrahlung 
emission, as both processes (p-p and free free emission) 
are proportional to the gas density squared (Katz & Wax- 
man 2008). As we show below, this implies a correlation 
between the radio luminosity and the X-ray luminosity, 
provided that the fraction of ICM thermal energy carried 
by CRs is roughly the same in different clusters. 

In our simple model the ICM is composed of a thermal 
gas of ionized Hydrogen, with temperature T and number 
density n, containing a population of proton CRs with 
a power law distribution e 2 dn/de = f3 COIC 3nT/2, and a 
magnetic field strong enough to make synchrotron emis- 
sion the dominant energy loss process for the radio emit- 
ting electrons. Our choice of the CR spectral index, for 
which the energy density per logarithmic energy interval 
is independent of energy, is consistent with the observed 
spectral index of the radio emission, L v oc v^ 1 (see also 
a detailed discussion of the Coma cluster observations 
in Kushnir & Waxman 2009b). The X-ray emission is 
produced in this model by thermal bremsstrahlung and 
has an emissivity (assuming the thermal Gaunt factor to 
be 1 and neglecting heavier than Helium elements, which 
can increase the bremsstrahlung emissivity by a few tens 
of percent): 

ex « fxixWraecV^T^n 2 , (2) 

2 A somewhat steeper correlation, Pi. 4 oc 1 2 4] 451 was "ob- 

tained by Brunetti et al. (2007). This is shown in magenta in 
figure 1. It is not clear to us how this fit was obtained. 
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where cr? is the Thomson cross section, a e is the fine 
structure constant and fx is a dimensionless factor of or- 
der unity which depends on the Hydrogen mass fraction 
X- The radio flux is produced by synchrotron emission 
of secondaries (electrons and positrons), which are the 
products of p-p collisions between the proton CRs and 
the thermal gas. We assume that the distribution of the 
radio emitting secondaries is in a steady state, where in 
the relevant energy bands the secondaries that are gener- 
ated lose all their energy to synchrotron radiation. This 
implies that the synchrotron emissivity per logarithmic 
frequency interval is 



1 



(3) 



Here ee e e is the energy production rate of secondaries 
per logarithmic secondary energy interval, given by 



eel 



o.i • /, 



pp 



(X)e 2 f nca™ 1 



pp 



(4) 



CTpp Cl w 40 mb is the cross section for inelastic collision 

and /p P (x) is a dimensionless factor of order unity, which 
depends on the Hydrogen mass fraction. 

Using equations (2), (3) and (4) we see that in this 
model the radio luminosity per logarithmic frequency in- 
terval is proportional to the bolometric X-ray Luminos- 
ity, 
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J 0.075a: 



T 



1/2 



1/2 _incl r 

PP JPP 



(x) 



CTT fx(x) 



~ 1.1 • lO-^core,-^ 1 ^, (5) 

where /3 CO ro,-4 = /^core/10 -4 , T\ = T/lOkeV and we used 
X = 0.75. If ,3corc does not vary considerably between 
clusters we expect a linear correlation between vL s J nc 
and T^ 2 L X . 

Figure 2 compares the prediction of eq. (5) with the 
distribution of the 17 clusters of Cassano et al. (2006) in 
the ln(j/L s / nc )-ln(T 1 / 2 L x ) plane (the values for T and 
Lx, as well as the measurement errors in T and Lx, are 
taken from Cassano et al. 2006). The error bars of the 
radio luminosity describe our estimate of the intrinsic 
scatter, AlnP = ln(1.7). A linear fit for the correla- 
tion between ln(j/L^ ync ) and ln(T 1 / 2 Lx) gives a slope of 
1.3 ± 0.7, consistent with the predicted slope of 1. The 
best fit obtained for a slope of 1 gives /3 CO rc,-4 — 2 (see 
eq. (5)). The fact that there is a scatter of about a factor 
2 around the correlation probably implies that there is a 
similar scatter in the value of /3 C oro among clusters. This 
is consistent with the scatter expected if the source of 
cluster CRs is acceleration in accretion shocks (Kushnir 
& Waxman 2009b, see discussion in § 6). 

Note that equation (5) can be applied to the ratio of 
the radio to X-ray surface brightness at any given point. 
For example, we give in table 1 the central (bolomet- 
ric) X-ray surface brightness, Fx, and the central radio 
surface brightness (taken as the maximal halo surface 
brightness), vFj*, ync , for the three clusters in Govoni et al. 
(2001), for which such data are available. The value of 
/3 corc derived from the data using equation (5) is consis- 
tent with ~ 2 • 10~ 4 . 

In order to compare the predicted correlation, eq. (5), 
with the correlation derived from the data shown in fig- 
ure 1, eq. (1), we need to relate the bolometric X-ray 
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Fig. 2. — The distribution of the 17 clusters of Cassano et al. 
(2006) in the \n(uL s J )-ln(T 1 / 2 L x ) plane. The error bars of 
the radio luminosity describe our estimate of the intrinsic scat- 
ter, AlnP = ln(1.7). A linear fit for the correlation gives a slope 
of 1.3 ± 0.7 (best fit shown by the solid line), and the best fit for a 
slope of 1 (dashed line) gives /3 CO r C ,-4 — 2 (see eq. 5). 



luminosity Lx to T and to Lx[o.i,2.4]- Lx and L X [o.i,2.4] 
are related by 



Lx ~ L X [0. 1,2.4] 



[ Xm " e-*/ 2 K (x/2)dx' 



(6) 



where x = hv/T and Kq is the zeroth order modified 
Bessel function of the second kind (we have included here 
the Gaunt factor corrections, which cannot be ignored at 
this level of approximation, see Drummond 1961). For 
the relevant cluster temperatures we may approximate 

L x ~ 3T®- 6 L X [o.i,2A]- (7) 

Next, it is well known that Lx and T are strongly cor- 
related (e.g. Markevitch 1998; Arnaud & Evrard 1999; 
Reiprich & Bohringer 2002), 

Lx ~ LxoTi L , (8) 

with oll in the range 2.5 — 3 (The deviation of the lumi- 
nosities from this phenomenological relation is remark- 
ably small, of the order of tens of percents). Using eq. (7) 
and eq. (8), eq. (5) gives 



1+; 



J X[0. 1,2.4] ' (9) 

Note, that the predicted slope of the correlation, 
dln(^ ync )/dln(L X [o.i,2.4]) = 1 + depends 
weakly on the value of at (within the range of 2.5 to 
3). Adopting Lxo = 3 • 10 45 and ul = 2.5, our predicted 
correlation, eq. (5), may finally be written as 

-l 



Pi 



1.4 



2.5 • 10 31 ixfo.i,2.4],45^corc,-4 erg S 



Hz" 



(10) 



Comparing eq. (10) with eq. (1) we find that the pre- 
dicted slope of the correlation is consistent with the ob- 
served slope, and that the normalization of the corre- 
lation is consistent with the observed one for /3 corc ~ 
2 • 10~ 4 . 
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Table 1. /3 ore derived from surface brightness data 



Cluster name 


F x a ' b 






/?corc,-4^ 




[erg cm -2 s -1 arcmin~ 2 ] 


[keV] [ 


erg cm -2 s — 1 arcmin - 2 ] 




A520 


1.6 • icr 12 


8.3 


9 ■ 10~ 17 


5.7 


A773 


3.5 -10- 12 


8.6 


5 ■ 10~ 17 


1.4 


A2744 


1.9 • 10~ 12 


10.1 


5 ■ 10~ 17 


2.4 



a Data taken from Govoni et al. (2001). 

b The central (bolometric) X-ray surface brightness. 

c The central radio surface brightness (taken as the maximal halo surface brightness). 
d Dcrived by using equation (5). 



Finally, let us determine the minimum magnetic field 
strength and the minimum time for variations of the mag- 
netic field and of the electron injection, for which our 
model assumptions of fast electron cooling dominated 
by synchrotron losses are valid. The two main energy 
loss mechanisms of the electrons are synchrotron emis- 
sion and inverse-Compton scattering of CMB photons. 
In order to ensure that the dominant process through 
which the electrons lose their energy is synchrotron ra- 
diation, the energy density in the magnetic field should 
be larger than the CMB energy density. This implies 

B > B CMB = (cWT^j) 1 ^ ~ 3.2(1 + z f M G. (11) 

This implies that the magnetic field energy caries a frac- 
tion 

e B ~ 1.7 • lQ- 2 T^nZ\{B/B CMB ) 2 (12) 

of the ICM thermal energy, where n_3 = n/10~ 3 cm~ 3 . 
€b values of order a few percent are reasonable, and ob- 
served in large-scale systems (Kim et al. 1989). Note, 
that as long as eg C 1 the magnetic field docs not affect 
the dynamics of the ICM. 

The lower limit on B sets an upper limit to the electron 
cooling time. The frequency at which an electron emits 
most of its synchrotron power is given by v = voj 2 , where 
vq = 3eB / (4irm e c) and 7 is the Lorentz factor of the 
electron. Thus, 

^- 10 ^iiiai) 1/2 (^)" /2 ' ( 13 ) 

and the cooling time of the electrons is 

/ X -1/2 / N -3/2 

tcool « 0.23 (^fe) (^) Gyr. (14) 

For our model to be valid, the time scales for significant 
changes in the secondary injection and in the magnetic 
field must be larger than - 0.25(B/B C mb) _3/2 Gyr. An 
estimate for these time scales is given in § 4, and it shown 
that they are much longer than the cooling time. Note, 
that since the energy of the secondary electrons is ~ 0.1 
of the energy of the parent CR protons leading to their 
production, and since the energy of CR protons in clus- 
ters is expected to reach <~ 10 18 eV (for details, see Kush- 
nir & Waxman 2009b), we expect the electron energy 
distribution to extend well above 7 ~ 10 4 . 

4. MAGNETIC FIELD EVOLUTION 

We next address the bimodality of the cluster radio lu- 
minosity distribution. Note, that since the cooling time 



of the CR protons at the relevant energies is very long 
compared to the Hubble time and their escape probabil- 
ity from the cluster is very low (for details, see Kushnir 
& Waxman 2009b), once CR protons are present within 
the cluster core in sufficient quantities to produce a RH, 
they remain there. Hence, the absence of RHs in some 
clusters must be due, in our model, to lower values of the 
magnetic field in their ICM. Observations suggest that 
clusters that contain RHs show merger activity (Ven- 
turi et al. 2007). This suggests the following simple sce- 
nario: ICM magnetic fields are enhanced by mergers to 
B 3> B CMB , and later decay to lower values, B < 1 /j,G. 
Decay to B < 1 = -Bcmb/3 implies a suppression 
of synchrotron luminosity by a factor of 10, rendering 
it undetectable. As explained in § 3, the magnetic field 
required to explain the radio luminosity, B > Bqmb, 
carries only a small fraction of the ICM thermal energy 
[cq. (12)]. Numerical simulations support the possibility 
of amplification of seed fields to such values by turbu- 
lence, generated in cluster mergers (e.g. Ryu et al. 2008). 
We expect that in this scenario the magnetic fields will 
decay on time scales similar to the few Gyr time scales 
in which turbulence decays (see e.g. Subramanian et al. 
2006). 

Let us consider next the magnetic field decay time, 
which is required to account for the bimodality of the 
radio luminosity distribution. Before analyzing quanti- 
tatively the observed bimodality, two points should be 
stressed. First, only clusters from a complete sample 
can be used for a statistical study. We therefore restrict 
our analysis only to the complete sample presented in 
(Brunctti et al. 2007). Note, that upper limits to the 
radio luminosity for clusters without RHs were obtained 
only for 20 of the clusters. Second, it is sensible to discus 
the bimodality in the Lx — Pi. 4 plane only where the ra- 
dio fluxes implied by the correlation are far above the de- 
tection limit (say a factor of 10). We thus restrict the dis- 
cussion to X-ray luminosities Lx[o.i,2.4] ^ 10 45 ergs _1 . 

In this restricted, complete sample there are 6 clus- 
ters for which the correlation holds, 1 cluster with a RH 
flux lower than implied by the correlation, and at least 
9 clusters with radio luminosity that is smaller than the 
luminosity implied by the correlation by a factor of > 10. 
These small numbers allow only a rough estimate of the 
decay time idec of the magnetic field: Since the life time 
of massive clusters is tnf c ~ 5 Gyr and only ~ 1/3 of 
the clusters have RHs, the typical decay time should be 
roughly t dcc ~ tiif e /3 ~ 1.5 Gyr. 

A note is in order regarding the claim of Brunetti et al. 
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(2007, 2008), that the "empty" region in the L x - Pi a 
plain (between the upper limits on Pi a for clusters with 
no RHs and the line defined by the Lx — Pi a corre- 
lation) implies that the magnetic fields must decay on 
0.1 — 0.2 Gyr timescales. First, we note that this state- 
ment was not supported by a quantitative statistical 
analysis. Moreover, the sample presented in (Brunetti 
et al. 2007) seems consistent with much longer decay 
times, of the order of few Gyr. Given the low statis- 
tics, the fact that there is at least 1 cluster with detected 
radio emission, which is smaller than implied by the cor- 
relation, is consistent with a scenario in which clusters 
spend similar times satisfying the correlation and pro- 
ducing radio flux lower than implied by the correlation 
(say 3 Gyr and 1.5 Gyr respectively), with tdcc ~ 1 Gyr. 

According to the scenario outlined above, the magnetic 
field is enhanced during mergers and later decays. In 
addition, an increase by a factor of order unity of the 
secondary injection rate (and X-ray luminosity) may take 
place during and following a merger, due to the increased 
densities and temperatures. These variations occur on 
time scales similar to, or larger than, the dynamical time 
scale of the mergers, idyn ~ 1 Gyr, which is much larger 
than the electron cooling time, given in cq. (14). This 
implies that our model assumption, that the electrons 
cool down on a time scale short compared to timescales 
over which the magnetic field or the CR electron injection 
change considerably, is valid. 

It should be emphasized here that the possible increase 
in secondary injection rate following cluster mergers is 
not expected to modify the ratio of radio to X-ray lumi- 
nosity given in eq. (5), and is not expected therefore to 
introduce a scatter to the predicted correlation, eq. (10). 
This is due to the fact that fast cooling of the secon- 
daries, t CO oi <C tdyn, ensures that they are in a quasi 
steady state for which eq. (5) holds. Note that the ratio 
given in eq. (5) may be modified if significant particle ac- 
celeration takes place in merger shocks. However, such 
acceleration is not expected to significantly modify the 
existing CR population, due to the low Mach numbers of 
merger shocks (see e.g. Ryu et al. 2003; Gabici & Blasi 
2003; McCarthy et al. 2007; Skillman et al. 2008; Kushnir 
& Waxman 2009b). 

5. OTHER MODELS 

In the previous sections we presented a simple model 
for clusters with RHs in which the radiation is gener- 
ated by secondary electrons radiating in strong magnetic 
fields, B > £>cmb- In this section we consider alterna- 
tive models, where either B < Bqmb or the radiating 
electrons are not secondaries produced by p-p collisions, 
and argue that such alternative models are unlikely to 
reproduce the observed correlation of radio and X-ray 
luminosities. 

In models where B < Bqmb (e.g Dennison 1980; Blasi 
& Colafrancesco 1999; Brunetti et al. 2001; Petrosian 
2001; Cassano & Brunetti 2005; Cassano et al. 2007) the 
synchrotron emission at radio frequencies strongly de- 
pends on the magnetic field value: P ra dio B^ p+1 ^ 2 cx 
B a+1 , where p rts 3 and a rts 1 are the spectral indexes 
of the electron population and emitted radiation respec- 
tively, implying P cx B 2 . Any such model must thus as- 
sume a tight correlation between the magnetic field value 
B and the X-ray luminosity, with an allowed deviation of 



a few tens of percents (to explain the small scatter in the 
correlation between the X-ray and radio luminosities). 
It is challenging to find a physical mechanism generating 
such a tight correlation between B and Lx (see however 
Dolag 2006). It should be added that for B < B C mb 
the fraction of energy carried by the magnetic field is 
small, e B = B 2 /(8w)/(3/2nT) - 1.7 x lQ-^nZ\T^ l B_ & , 
where £>_6 = B/fiG, and thus large variations are not 
constrained by energetic considerations. 

Let us consider next B > £>cmb models in which 
the radiating electrons are not of secondary origin. For 
B > Bqmb , the radio luminosity and the rate of produc- 
tion of radiating electrons are related by ee e E rts 2ve s J nc 
(assuming that the time scale for variations in ee e e is 
longer than the cooling time). Since the synchrotron 
emission of secondary electrons accounts for the observed 
radio luminosity for f3 COIC ~ 10 , the generation rate of 
high energy electrons required to account for the radio 
emission is given by (see eq. 4) 

se e E ~ 10- 5 n 2 Tc4 n p el . (15) 

In a model where the radiating electrons are secondaries 
from p-p collisions, the scaling et% cx n 2 T is natural, and 
the normalization, eel/n 2 Tccr™° 1 w 10~ 5 corresponding 

to /3 CO rc ~ 10~ 4 , is naturally obtained for a reasonable 
value of the fraction of accretion shock energy converted 
to CRs (see discussion in § 6). It is challenging to ac- 
count for the scaling and normalization of the high energy 
electron production rate in models where the radiating 
electrons are not secondaries. 

Finally, we note that models, where the radiating elec- 
trons are not secondaries, require /3 CO rc <C 10~ 4 in order 
that secondary electrons would not dominate the radio 
flux. For example, this requirement is not satisfied in 
the model presented in (Brunetti et al. 2001; Cassano 
et al. 2007; Brunetti & Lazarian 2007), where the radiat- 
ing electrons are not secondaries but rather preexisting 
nonthermal CR electrons reaccelerated by cluster turbu- 
lence. In this model the CR density assumed in order 
to reproduce the observed radio luminosity (by emission 
from turbulence accelerated electrons) is a few percent, 
corresponding to /3 cor o ~ 10~ 2 . This implies that syn- 
chrotron emission from secondary electrons would dom- 
inate the radio luminosity in such models, and would 
exceed the observed radio luminosity by <~ 2 orders of 
magnitude. 

6. DISCUSSION 

We have presented (§ 3) a simple model that explains 
the observed radio emission from clusters as synchrotron 
radiation from secondary electrons that are produced by 
p-p interactions and lose their energy by emitting radi- 
ation in a strong magnetic field, B > {^oT^j^g) 1 / 2 ~ 
3 fiG. The radio luminosity is given in this model by the 
rate of production of high energy electrons (eqs. 3, 4). 
This, combined with the observed correlation between 
the X-ray luminosity and the temperature of clusters 
(eq. 8), naturally reproduces the observed radio and X- 
ray luminosity correlation, Pi. 4 cx L^J 1 2 4 j (see cq. 1), 

provided that the fraction of ICM thermal energy carried 
by CRs per logarithmic particle energy interval is roughly 
the same for all clusters, and equal to /3 C orc "~ 2-10- 4 (see 
eq. 10). We have shown (§ 5) that alternative models, 
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where either B < Pcmb or the radiating electrons are 
not secondaries produced by p-p collisions, arc unlikely 
to reproduce the observed tight correlation of radio and 
X-ray luminosities. 

The magnetic field implied by our model, B > Bcmb, 
carries at least a few percents of the thermal energy of 
the ICM, e B > 1CT 2 [Eq. (12)]. This is much higher 
than the fraction /3 C ore,tot. of thermal ICM energy carried 
by the CRs: for a e 2 dn/de oc e° CR spectrum we have 
/^core, tot. — 10/3 C oro ~ 10~ 3 < O.les. This is an indication 
that these two nonthermal components are governed by 
different processes (as we suggest below, the CR protons 
are generated at the accretion shocks, while the magnetic 
fields are amplified during mergers). 

Clusters that contain RHs show merger activity (Ven- 
turi et al. 2007). In order to explain the bimodality of 
clusters in the P1.4— £x[o.i,2.4] plane reported by Brunetti 
et al. (2007), we suggested in § 4 that the magnetic field 
in clusters is amplified following mergers to B ^> Bcmb, 
and later decays on time scales of 1 Gyr to lower values, 
B < 1 fj,G. Numerical simulations support the possibil- 
ity of amplification of seed fields to such values, which 
correspond to few percent of equipartition (eq. 12), by 
turbulence generated in cluster mergers (e.g. Ryu et al. 
2008). The time scale for magnetic field decay, implied 
by the Pi. 4 — Lx[o.i,2.4] distribution, is similar to the 
few Gyr time scale of turbulence decay (theoretical ar- 
guments supporting magentic field decay on a time scale 
similar to that of turbulence decay may be found, e.g., in 
Subramanian et al. 2006). We emphasis that the obser- 
vational fact that only ~ 1/3 of the clusters host a RH 
(and these clusters show merger activity) is an indication 
for deviation from the minimum energy state (for which 
the energy density in magnetic fields equals the CRs en- 
ergy density), for clusters that host a RH. Note, that for 
the other ~ 2/3 of the clusters the minimum energy state 
might hold, since in our scenario the CR protons' energy 
density is not affected significantly by mergers, but the 
energy density of the magnetic field decreases (after the 
merger) by a factor of ~ 10 (at least, to match the up- 
per limits), which makes the energy densities of the two 
nonthermal components roughly equal. 

The observed tight correlation between P 1A and L x 
suggests that the fraction of ICM thermal energy car- 
ried by CRs is indeed uniform among clusters and given 
by /3 C oro ~ 2 • 10~ 4 . Various sources of CR were con- 
sidered in the literature, including active galactic nuclei 
(e.g. Katz 1976; Fabian et al. 1976; Fujita et al. 2007; 
Blasi et al. 2007), dark matter bow shocks (e.g. Bykov 
et al. 2000), ram-pressure stripping of infalling galaxies 
(e.g. de Plaa et al. 2006), supernova (SN) explosions (e.g. 
Volk et al. 1996) and shock waves associated with the 
process of large scale structure formation (e.g. Loeb & 
Waxman 2000; Fujita et al. 2003; Berrington & Dermer 
2003; Gabici & Blasi 2003; Brunetti et al. 2004; Inoue 
et al. 2005). The resemblance between cluster accretion 
shocks and collisionless non-relativistic shocks in the in- 
terstellar medium (Loeb & Waxman 2000), which are 
known to accelerate a power-law distribution of high en- 
ergy particles, suggests that accretion shocks may be the 
source of cluster CRs (Note that acceleration in merger 
shocks is not expected to significantly modify the ex- 
isting CR population, due to the low Mach numbers of 



merger shocks, e.g. Ryu et al. 2003; Gabici & Blasi 2003; 
McCarthy et al. 2007; Skillman et al. 2008; Kushnir & 
Waxman 2009b). 

In order to account for the observed radio luminosity, 
the fraction i] p of accretion shock energy deposited in 
relativistic CR protons should be a few percents. This 
conclusion is reached by noting that /3 C oro.tot. is related 
to i] p by r\ v <~ 10/3 CO ro,tot. (Pfrommcr et al. 2007; Kush- 
nir & Waxman 2009b) . The latter relation is determined 
mainly by the reduction of the thermal energy fraction 
carried by CRs as the shocked gas is adiabatically com- 
pressed falling into the cluster's center. The compression 
of the gas leads to a reduction of the CR energy fraction, 
oc p^ 1 / 3 , where p is the gas density, due to the softer 
equation of state of the relativistic particles (see Kushnir 
& Waxman 2009b, for a detailed discussion). The com- 
pression also increases the CR energy by a factor ~ 10, 
such that T] p estimated above is relevant for 10 GeV CR 
(note, that the parent CR energy is <~ 10 times the en- 
ergy of the secondary produced by an inelastic nuclear 
collision). 

The accretion shock origin of cluster CRs is supported 
by both the inferred value of r\ p and the inferred scatter 
in /3 coro : The scatter in /? C ore,tot., in a model where CRs 
are accelerated in accretion shocks and advected towards 
the cluster center with the infalling gas, yields a scatter 
of roughly a factor of 2 in the value of /3 C orc among dif- 
ferent clusters (Kushnir & Waxman 2009b), consistent 
with the scatter implied by the Pi. 4 — Lx relation; The 
implied acceleration efficiency of protons, r\ p ~ 5 per- 
cent, is similar to the acceleration efficiency of electrons 
in the accretion shocks, implied by the nonthermal HXR 
emission detected in several clusters (Kushnir & Wax- 
man 2009b). 

We next discuss the possibility that the origin of ICM 
CRs is SN explosions (e.g. Volk et al. 1996). In this 
model, the observed Iron abundance of the ICM (e c ; ~ 
1/3 relative to the solar value [Be]©) is used to estimate 
the energy released by SN explosions, assuming some ra- 
tio between the energy released in SNe and the Iron mass 
ejected per explosion, (Esn /SMp e ). Assuming CRs are 
also injected in such explosions with acceleration effi- 
ciency per logarithmic energy bin of Z?^, the ratio of 
the energy density of such CRs to the thermal energy 
density of the ICM per logarithmic energy bin is given 

by 

E SN gSNi? ( -I c\ 

) ~5MF~ e CR ' (16) 



SN 



pm p 



;£ c i[Fe 



(3/2)T 

The last factor P coo i represents the energy loss of these 
CRs as they leave galaxies into the ICM. The ratio be- 
tween f3 s ^ e and the corresponding value for accretion 
shocks /S^ore i s given by 



SN 

core ^ ^ q f^CR p 



3acc 
J core 



JSN 
?&cc 



( Esn 



cool 



5M 



Fe 



Vio 5i cr g ; vo.im ; ^ ' ^ 

where is the acceleration efficiency per logarithmic 
energy bin in accretion shocks, and we included the re- 
duction of the thermal energy fraction carried by CRs as 
the shocked gas is adiabatically compressed falling into 
the cluster's center. The adopted value for {Esn /5Mp e ) 
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is probably an upper limit, since it is derived by assuming 
that all of the Iron is generated in type II SNe (see how- 
ever Hamuy 2003, for the large uncertainty in this value). 
However, some of this Iron can be produced in other pro- 
cesses, especially type la SNe with SMp e ~ 0.7M Q . Fur- 
thermore, -Fcooi is probably a small number, such that 
if the acceleration efficiency in both SNe and accretion 
shocks is similar, than we expect accretion shock CRs to 
dominant those produce in SNe. 

In our analysis we have so far neglected the possible 
effects of CR diffusion. If 100 GeV CR diffusion is sig- 
nificant over scales > 100 kpc, CRs may "escape" the 
infalling gas, possibly leading to f3 COIC ^ot./Vp ^ 1/20. 
A discussion of the transport of CRs in the collisionless 
cluster plasma, a subject which is poorly understood, is 
beyond the scope of this paper (note, that CR diffusion is 
poorly understood also in the Galaxy, where much more 
data are available on the CRs and on the magnetic field 
configuration, e.g. Ptuskin et al. 2006). We note, how- 
ever, that if diffusion were important, allowing most of 
the CRs accelerated in the accretion shock not to reach 
the cluster core, then j3 corc would have been likely to vary 
considerably between clusters. The fact that we observe 
a tight correlation between the radio and X-ray luminosi- 
ties indicates that this is not the case. 

According to our model, the clusters' radio luminos- 
ity is given by the energy production rate of charged, 
e , secondaries. This implies that the high energy 



gamma-ray luminosity produced at the cluster core 
by the decay of neutral, 7r° secondaries, is propor- 
tional to the radio luminosity and should be given 
by uLZ~ Tay ~ 3 • 10 41 ixj r . 12 . 4]i45 ergs- 1 . The predicted 
core luminosity will be difficult to detect with existing 
and planned instruments. Note, however, that high en- 
ergy inverse Compton gamma-ray emission, from elec- 
trons accelerated at the accretion shocks, should be de- 
tectable by FERMI and possibly by existing Cherenkov 
telescopes (see detailed discussion in Kushnir & Waxman 
2009b). 

Finally, a cautionary note should be made regard- 
ing predictions for low frequency radio observations 
of clusters (e.g. with ALMA, LOFAR). Under our 
model assumptions, the characteristic energies of elec- 
trons dominating the emission at low frequencies are ~ 

1.4(!/ radio /100MHz) 1 /2( J B/S C MB)" 1/2 GeV, which are 
the products of CRs that were accelerated to ~ 1 GeV 
energies. Since the spectrum of CRs is poorly understood 
at energies e ~ m p c 2 , the extrapolation of the radio spec- 
trum from the currently observed high frequencies to the 
relevant low frequencies is uncertain. 
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